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Abstract 

Conformational properties of star-shaped polymer aggregates that carry attractive end-groups, 
called telechelic star polymers, are investigated by simulation and analytical variational theory. 
We focus on the case of low telechelic star polymer functionalities, / < 5, a condition which 
allows aggregation of all attractive monomers on one site. We establish the functionality- and 
polymerization-number dependence of the transition temperature from the "star burst" to the 
"water melon" macroparticle structure. Extensions to telechelic stars featuring partially collapsed 
configurations are also discussed. 

PACS numbers: 61.25.Hq, 82.70.Uv, 36.20.Ey, 61.20.Ja 
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Self-organizing soft materials are relevant in developing novel macromolecular compounds 
with peculiar structural and dynamical properties, which are related to mesoscopic aggre- 
gation as well as intra- and inter-molecular association Progresses in the synthesis 
of chains with attractive end-groups, called functionalized or telechelic polymers, opened 
the way for the subsequent synthesis of telechelic star polymers with attractive polar end 
groups, telechelic associating polymers with hydrophobic terminal groups^ssociating poly- 
electrolytes in homogeneous solutions, and telechelic planar brushes 0,0,0,0. A telechelic 
star, which is the subject of the work at hand, consists of / chains with one of their ends 
chemically attached on a common center, whereas the attractive ends are free at the other 
end of the molecule. Block copolymer stars with a hydrophilic core and a hydrophobic corona 
0, ^ also bear similarities to telechelic stars, the two systems becoming identical when the 
length of the hydrophobic block reduces to a single monomer. The thermodynamics and the 
structure of planar telechelic brushes have been recently analyzed also from a theoretical 
point of view, shedding light into the quantitative characteristics of both the conformations 
and the interactions of the same . Further theoretical approaches have been devel- 
oped to describe fiower-like micelles with hydrophobic terminal groups that self-assemble in 
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water. Such aggregates show a characteristic 'bridging attraction' 
to a liquid- vapor phase transition [13] . The interesting feature of telechelic star polymers 
is the possibility of attachment between the attractive terminal groups. In dilute solutions, 
this gives rise to intra-molecular association as well as inter-association between micelles, 
depending on the details of the molecular structure. 

Examples of low-functionality telechelics, very similar to the ones considered in this work, 
are mono-, di-, and tri-cj-zwitterionic, three-arm star symmetric polybutadienes. Experi- 
ments have shown that these self-assemble into distinct supramolecular structures, includ- 
ing collapsed, soft-sphere conformations ^ 14|. In particular, using low-angle laser light 



scattering and dynamic light scattering, it was found that samples with three zwitterion 
end groups present a low degree of inter-association between macromoleculeSjShowing in- 
stead a preference for intra-association and formation of collapsed soft spheres 0, 3]. X-ray 
scattering and rheological experiments support the conclusion that this tendency persists at 
higher concentrations, all the way into the melt ^. There, the formation of transient gels 
has been found for the case of two- and three-zwitterion macromolecules, with the network 
*a.aote.stic. depending on .Cecula. we.g.. of tKe „n. Q. 
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The conformations of telechelic micelles are evidently determined by the competition 
between entropic and energetic contributions. A detailed investigation, by theory and simu- 
lation, of the mechanisms leading to the formation of collapsed soft spheres is, however, still 
lacking. In this Letter we perform extensive computer simulations, accompanied by a scaling 
analysis of the free energies of candidate structures of telechelic micelles with small function- 
ality, / < 5. We find that at high temperatures the system assumes the usual star burst (sb) 
configuration. On the contrary, at sufficiently low temperatures, the end-monomers attach 
to each other and the micelles assume an overall closed configuration, akin to the collapsed 
soft spheres conjectured in the experimental study of Ref. j^. Due to the peculiar shape 
of these aggregates, featuring two points of aggregation, one at each end, we term them 
"water melons" {wm). For intermediate temperatures, the macromolecules exhibit several 
configurations that correspond to a partial assembling of the terminal groups. The number 
of chains connected at their ends depends on the temperature T, the functionality / and the 
degree of polymerization N of the chains. In particular, at a given, sufficiently low T, the 
wm-configuration becomes more stable with respect to the sb-one with decreasing and 
increasing /, as will be demonstrated in what follows. Our theoretical predictions are found 
to agree well with simulation results. 

We employed monomer-resolved molecular dynamics (MD) simulations to examine the 
conformations of isolated telechelic micelles. The monomers were modeled as soft spheres 
interacting by means of a truncated and shifted Lennard- Jones potential, ^j(r) (length 
scale (Tlj, energy scale e), the truncation point being at the minimum, = 2^/^(Tlj, 
rendering the interaction purely repulsive. The end-monomers, on the other hand, interact 
with each other by means of the full Lennard- Jones potential, Vi^j{r) that has a minimum 
value VLj(rmin) = The chain connectivity was modeled by employing the finite extension 
nonlinear elastic (FENE) potential Q, using the same parameter values as in the simulation 
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of normal star polymers |17l . Il8| . The monomer mass m, the energy e and the length ulj 
are assigned the value unity, defining thereby the characteristic time r = ma^^j/e and the 
dimensionless temperature T* = k-QT/e, where fee is Boltzmann's constant. The equations 
of motion were integrated using the velocity form of Verlet's algorithm 2, whereas the 
temperature was fixed by applying a Langevin thermostat 20]. The timestep was At = 
10~^r, with a total of 2 x 10^ timesteps used for equilibration and 10 independent runs of 
5 X 10^ timesteps each to gather statistics. The characteristic quantities measured were the 
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FIG. 1: The probability distribution function P{Et) obtained from simulations of / = 3 telechelic 
micelles carrying N = 10, 50 monomers per arm. Note the logarithmic scale. 



radius of gyration Rg of the molecules, the radial distribution function gt{r) between the 
terminal monomers and the expectation value Et of the interaction between the terminal 
groups, with Et = {J2i=i X]j>j ^jd^j ^ where (■ ■ ■ ) denotes a statistical average and 
tj stands for the position vector of the end-monomer of the i-th chain. 

We considered temperatures between T* = 0.01 and T* = 1.2. For T* > 0.2, we found 
that the micelle configuration is similar to that of a star polymer without attractive ends, i.e., 
the s6-configuration. A morphological transition of the micelles takes place upon lowering 
the temperature. In Fig. Q we show the probability distribution function P{Et) for the 
end-monomer interaction energy Et for the case / = 3, for different T*- and A^-values. 
Let us focus on = 10. At T* = 1.0, Fig. Ufa), P{Et) shows almost all the events 
around Et = 0, meaning that in the vast majority of configurations the end-monomers are 
far apart, so that the attractions are vanishing. For T* = 0.2, Fig. Wljo), P{Et) takes a 
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bimodal form with two peaks, one at Et = and one at Et = —e. The latter corresponds 
to a conformation in which two chains are end-attached, with their terminal monomers 
at a distance r^^^^, whereas the third chain is still free (intermediate configuration). Note 
that the molecule does fluctuate between the two conformations, as witnessed both by 
the bimodal character of the probability distribution and by the non-vanishing values of 
P{Et) for —e < Et < 0. Upon further lowering of the temperature at T* = 0.1, Fig. 
life), P{Et) develops a peak at Et = —3e, corresponding to a configuration in which all 
three terminal monomers are confined at a distance r^[^ from each other. This is the state 
in which the micelle has aggregation points at both ends, which we termed 'water melon' 
(?i'm)-conformation. The trend described above persists for = 50. However, in each panel 
we can notice a higher probability to have association among chains by decreasing A^. We 
further found a considerable sharpening of the peak at Et = —3e for lowered temperatures, 
since fluctuations around the ground state are getting more suppressed. 

Similar results have been found for the cases / = 2, 4, and 5. At fixed A^, the probability 
of chain association grows increasing /. The minimum value of Et for / = 2 is —e and 
for / = 4 it is —6e, corresponding to a wm-configuration in which every end-monomer is 
confined at a distance r^[^ from every other one, i.e., an arrangement at the vertices of a 
regular tetrahedron of edge length rmin. Due to geometrical constraints, for / = 5 it is not 
possible to have all five end-monomers at contact with each other; the minimum value for 
Et is —9e in this case, corresponding to the arrangement of the terminal segments on the 
vertices of two regular tetrahedra of edge length rmin, which share a common face. 

The temperature dependence of the gyration radius Rg for fixed A^ = 10 and different 
/-values is shown in the main plot of Fig. |21 At high temperatures, Rg has a plateau 
that corresponds to s&-case and scales as f^/^N^/^ 21 1. As T* is lowered, we find a rapid 
decrease of Rg within a narrow temperature range and a saturation to a lower plateau 
value that corresponds to the size of the wm-configuration. The attachment of the terminal 
monomers leads to a shrinking of the molecule, in agreement with the experimental findings 
of collapsed soft-sphere conformations 0| . The wm-configuration persists to higher T*-values 
upon increasing the functionality /, a trend that can be attributed to the increasingly strong 
attractive-energy contributions to the W7m-free energy as / grows. In the inset of Fig. |21 we 
show the terminal-segment radial distribution function for the various /-values at T* = 0.1, 
where the wm-conformation is stable and compare it with the one obtained for / = 3 
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FIG. 2: (Color online) Temperature dependence of the gyration radius Rg obtained by MD of 
telechelic micelles with = 10. From top to bottom: f = 5 (black), f = ^ (blue), / = 3 (green), 
and / = 2 (red). Inset: the end-monomer radial distribution function gt{r) obtained from MD at 
T* = 0.1 (solid curves); the sequence from top to bottom and the color coding are the same as 
those in the main plot. The dashed line is gtir) for / = 3 and T* = 1.0 for comparison. For clarity, 
each curve has been multiplied by a different constant. 



at T* = 1.0, where the micelles assume a s6-conformation. The aggregation of the end- 
monomers is clearly witnessed by the high peak at r = rmin, which is present at T* = 0.1, 
as opposed by the fiat shape of gt{r) at T* = 1.0. Whereas for / = 2, 3 and 4 a single 
accumulation peak can be seen, a double peak is present for the case / = 5. This feature 
arises from the arrangement of the five terminal monomers at the vertices of two tetrahedra 
that share a common face. The most distant vertices of the two are separated by a distance 
■\/8/3rinin, at which indeed the second peak in gt{r) shows up. 
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For the theoretical analysis, we consider the sb- and w;m-configurations and put forward 
a scaling theory, to determine the free energy F of each and find the stable one (lowest F- 
value). The excluded volume parameter v is set to unity, consistent with the choice ctlj = 1 
above {v = cr'lj)- In the case of the sfe-configuration, the free energy is entirely entropic 
in nature and can be expressed as a sum of the elastic and excluded-volume contributions. 
For a single chain, minimization of this sum with respect to the chain radius Rg leads to 
the scaling laws Rg ~ N^/^ and F ~ k^TN^^^. For high-functionality star polymers, Flory 
theory can be improved in order to better take into account the inter-chain correlations by 



invoking the blob model of Daoud and Cotton |21l . |22| . For the small /-values at hand, 
we envision instead each chain as a blob of radius ~ Rg that occupies a section of space 
delimited by a solid angle 47r//. Accordingly, the total free energy of the 56-configuration is 
approximated by Fsb = k-^T fN^/^ . The attractive contribution of the terminal monomers 
is very small, since the latter are far apart in the sfe-state, and it can thus be ignored. 

In the t/;m-configuration, there are two accumulation points: the first is the point in 
which the chains are chemically linked and the second is at the other end, at which the end- 
monomers stick together. In order to properly take into account the monomer correlations 
within the macromolecule ^Sf , we emplgv a blob model to estimate the excluded-volume 
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contributions to the free energy Fwm |lfll.l24|. A schematic representation of the water melon 
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and the associated blob model are shown in Fig. El the wm is thereby modeled as a double 
cone. With tq and I being the radius and height of the cone that contains Ub blobs of a 
single chain, we obtain 

TP q -I- /2 
"""^ :/^ + 2M + E,,,,. (1) 



The first term at the right-hand side of Eq. above is the stretching contribution of / 



chains, each having an extension \/rQ + P. The second is the excluded-volume cost and 
arises from the total number of 2fnb blobs, each contributing an amount k^T [l^ Q ■ 
Finally the third term is the attractive energy contribution arising from the total number 
of contacts between terminal monomers. Accordingly, -Eattr = ~f{f ~ 1)/(2T*) for / = 2, 
3, and 4, whereas i^attr = — 9/T* for / = 5. The blobs are close-packed within each cone, 
hence: 

' „-f.V l + sinW2) y 

co8(a/2) Z.l^i_si„(a/2)J • ^ > 

where a = arctan(ro//) is the cone opening angle and D is the diameter of the first, smallest 
blob, taken equal to the monomer size. The polymerization is obtained by summing over 
the monomers of all blobs of a chain, taking into account that a blob of size b contains 
Nb ~ b^^^ monomers. This yields 

"6-1 /. I . / /r,\ \ (5n/3) 



"6 — / -, 
n=0 ^ 



sin(a/2) 



(3) 



We solved Eqs. (j2I) and Q numerically to express Ub and / in terms of tq and then we 
minimized numerically F^jm with respect to tq using Eq. (^. We found that the cone angle 
a between two arms increases rapidly with for small A^-values: for 5 < A^ < 20 it changes 
from 20 to 40 degrees. This trend reflects the physical difficulty for the chains to form 
an aggregate for short chains, due to strong steric hindrance at the star center. Once the 
excluded volume interactions are balanced, the angle increases slowly, by about 5 degrees 
for 25 < A^ < 400. 

Comparing the minimized free energies -F,„rp, with Fgb, we obtained the transition temper- 



ature between sb and wm as a function of A^ 



26|. In Fig. m we show the T vs. A^ transition 



lines, separating the wm-state (below) from the sp-one (above). The 'critical number' Nc 
increases with decreasing T* and the transition temperature increases with /. Stated oth- 
erwise, for fixed A^ the theoretical model gives evidence to a stronger stability of the wm 
configuration on increasing / and lowering the temperature. On the other hand, keeping / 

8 




25 50 75 100 125 150 175 200 225 250 275 300 

N 



FIG. 4: Temperature versus N diagram corresponding to the sb- and wm-morphological transitions 
of telechelic micelles. Main figure: theoretical results. Inset: transition lines in simulations, 
determined via the T* dependence of Rg . From bottom to top (full/empty symbols): / = 2 
(circle) , / = 3 (square) , / = 4 (triangle up) , / = 5 (triangle down) . 



and T* fixed, the s6-configuration is stable above a certain value of Nc. The dependence 
of the repulsive contribution in the free energy on the number of monomers is not triv- 
ial: increasing A^, the number of blobs rii, increases and so does the contribution due to 
the excluded volume interaction. The entropic contribution is proportional to Tq, which 
increases with A^, and to The balance between all these terms in Eq. determines 

the stable micelle conformation: it is physically more expensive for the molecule to assume 
a wm-configuration the lower / and the higher is. Nc decreases rapidly with T* and the 
wm-state becomes unstable at T* > 0.3 for all /-values considered. In the inset of Fig. I^we 
show the corresponding MD-results: notice the nice agreement with the theoretical predic- 
tion trends. The value of transition temperature is overestimated in theory, reflecting the 
mean-field nature of the latter. There exist also configurations intermediate to the wm- and 
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s6-ones, in which only a partial association of arms takes place, akin to those seen for block- 
copolymer stars with / = 12 in the simulation study of Ref. Such configurations are 
also expected to be dominant as / further increases, due to the very high entropic penalty 
inherent to a complete intramolecular association. These conformations can be analyzed 
along the lines put forward in this work, since intermediate states can be looked upon as a 
combination of smaller water melons (work in progress). 

We have studied the dependence of the equilibrium conformations of telechelic micelles 
on the temperature, arm number and functionality of the same. The system presents a 
stable, water-melon configuration for low temperatures that opens up as T grows. The 
dynamical properties of such aggregates in solution will have unusual characteristics: at low 
temperatures and near the overlap concentration, various types of gels can appear, owing 
their viscosity either to transient network formation or to the existence of entangled loops. 
An understanding of the changes in the micelle conformations is a relevant starting point, 
aiming at gaining control over the supramolecular structural and dynamical properties of 
the system [2^, and targeting specific apphcations. 
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for helpful discussions. CM. thanks the Diisseldorf Entrepeneurs Foundation for financial 
support. This work was supported by the DFG within the SFB-TR6 and by the Marie Curie 
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